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prove analytically that the maximal classical achievable fidelity (which is 1/2 
without squeezing, for s = 1) is given by y/s/(l + s), vanishing when the degree 
of squeezing diverges. For mixed states, as well as for general distributions of 
displacements, we reduce the determination of the benchmarks to the solution 
of a finite-dimensional semidefinite program, which yields accurate, certifiable 
bounds thanks to a rigorous analysis of the truncation error. This approach may 
be easily adapted to more general ensembles of input states. 
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1. Introduction 

The storage and entanglement assisted teleportation of quantum states are two 
of the central primitives of Quantum Information Science. They have by now 
been accomplished with increasing precision in various experimental settings, with 
outstanding examples in the 'continuous variable' regime, adopting light modes 
or collective spins of atomic ensembles to, respectively, carry and store quantum 
information [B HI HI IS], or relying on motional atomic degrees of freedom [5J [7]- 
The need to certify success in such experiments, and to justify the use of the term 
quantum in setups such as 'quantum memories' and 'quantum teleportations', requires 
theoretical benchmarks which bound the performance of purely classical schemes 
[H [lOl [Til H2]- Here, "classical" refers to protocols where the quantum system 
is measured and later re-prepared from information obtained in the measurement, 
which is in turn stored or transmitted by classical means. The fact that there are 
limitations to such 'measure-and-prepare' schemes immediately follows from the no- 
cloning principle. A precise notion of these limitations, however, depends on the figure 
of merit (usually the fidelity) as well as on the prior distribution, i.e., on the ensemble 
of quantum states to be stored or transmitted. 

If, for a (i-dimensional Hilbert space, the input ensemble is comprised of all 
the pure states distributed according to the Haar measure, the average fidelity 
achievable by a classical scheme is 2/[d + 1), dropping to zero with increasing 
dimension d [131 114). Clearly, for continuous variable systems (where d = oo) not 
all pure states are experimentally accessible (as they constitute a set with infinitely 
many real parameters and, in principle, unbounded energy). The typical continuous 
variable implementations via electro-magnetic field modes or collective fluctuations 
in atomic ensembles favor Gaussian states of which coherent states are the simplest 
representatives. It was proven in [5] that the optimal fidelity for classical schemes is 
1/2 if the input ensemble is made up of coherent states taken from a flat distribution in 
phase space. This shows that such restrictions on the 'alphabet' of input states allows 
for classical schemes to achieve finite average fidelities, even though the dimension of 
the Hilbert space is infinite. 

The present work deals with input ensembles including Gaussian squeezed states, 
for which one expects, in general, the constraint on classical schemes to become more 
and more severe as the degree of squeezing is increased. This study finds its motivation 
in recent and ongoing experimental attempts to teleport and store squeezed quantum 
states [15j [HI [17l [18] : we will provide a means for certifying success in experiments. 

Although the tools we shall develop are suitable for more general applications, 
encompassing non-Gaussian states or even finite dimensional systems, the focus of the 
present paper will be Gaussian input ensembles with a fixed degree of squeezing s, 
(s being the factor by which the variance of one of the two canonical quadratures is 
reduced with respect to the vacuum level). This complements the results of Ref. |llj . 
where s is assumed to be unknown, and generalizes the results of Ref. [9] on coherent 
states (special case s = 1 in our notation). In brief, we have been able to obtain the 
following results: 

(i) For an ensemble of squeezed coherent states where, apart from the degree of 
squeezing, no a priori information is given about orientation and displacement, 
the maximal classical achievable fidelity is yfsj (1 + s). The optimal measure-and- 
prepare scheme achieving this bound is realized by heterodyne detection followed 
by preparation of a coherent state. 
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(ii) For an ensemble of squeezed states as in Q, but where each state is subject to 
additive Gaussian noise, analytical bounds are derived for the fidelity as well as 
for the maximal overlap achievable by classical schemes. 

(iii) For ensembles with a random orientation (phase covariance) but arbitrary 
distribution of displacements the maximal overlap is shown to be computable 
by means of a semidefinite program (SDP). Together with a rigorous bound on 
the truncation error in Fock space this yields reliable benchmarks. 

(iv) The SDP is applied to various cases of pure state ensembles (with and without 
Gaussian displacement distribution) and mixed state ensembles for parameter 
regimes relevant to present and future experiments. 

Further technical details can be found in the Appendices, where we present a 
characterization of (time-reversible) covariant channels and a discussion of the Choi- 
matrix formalism in infinite dimensions. 



2. Notation and conventions 

In this Section we briefly recapitulate some basic notation and useful concepts. For a 
more detailed exposition the reader is referred to Refs. [l!?l 1201 12T] . 

Throughout the paper, we will deal with bosonic systems of n modes, each 
of which is assigned to an infinite-dimensional Hilbert space TL — -^(R). 58 (W) 
and €i(H) will denote the set of bounded operators and the set of trace-class 
operators, respectively. Let us arrange the In canonical operators in a vector 
R = (Xx,Pi, . . . , X n , P„) T such that [Rj,R k ] = io^l, o = i ©™ =1 cr y being the 
anti-symmetric symplectic form (while o~ y is the two-dimensional y-Pauli matrix). 

As is customary, the Weyl displacement operator W% will be defined as W$ = 
e i ^ crR , for f £ E 2 ' 1 so that W^W V = W v W^- l ^ m K The 'characteristic function' 
Xa(£) of an operator A e is defined as 

Xa(£.) — Tr [W{A] . (1) 

In turn, the operator A is determined by its characteristic function according to the 
Fourier- Weyl relation 

which leads to the useful Parseval relation for A\,Ai G £i(TC® n ): 



Tr 



Am] = (^) n J n2n xa 1 (0xa 2 (£)d 2 "e . (3) 

For a density operator p, we define a 'covariance matrix' (CM) j p with entries 

(l P ) 3 k = Tr [{Rj,R k }p} - 2Tr [R jP ] Tr [R k p] 

and a vector of first moments d p with entries (d p )j = Tr [Rjp\. A state p is said to be 
Gaussian if its characteristic function \p is a Gaussian, reading 



x P (0 



e 



i£-crd p -(cr£)-7p(<7?)/4 



The vacuum |0) (0| is a Gaussian state with 7|o)(o| = 1 an d ^|o)(o| — 0- 
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3. Setup and figures of merit 

Our goal is to quantify the limitations for measure-and-prepare schemes T, i.e., 
elements of the set £ 3 T of entanglement-breaking channels |22[ 123] . when acting 
on an ensemble of input states characterized by a set of parameters w. Consider 
for instance an ensemble of pure squeezed Gaussian states with CM 

_ / cos 9 sin 9 \ / s \ / cos# — sin# \ 
7= ^ -sm6> cosfl ) \ l/s ) \ sin6> cos(9 )' ' 

and displacement £ and thus lo = (s,£,0). In order to fix a useful figure of merit 
we have to choose a functional F{T{p UJ ), p u ) which (i) measures, in some sense to 
be specified, the ability of T to 'preserve' the state p u when applied to it and (ii) 
can be determined in experiments. Based on this choice, we can then either quantify 
the worst case performance or the average case performance of a channel T, where 
the latter depends on an a priori distribution q(uj) over the parameter space. The 
corresponding benchmarks are then obtained by taking the supremum over all T G £ 
leading to the definitions: 

F (T) := inf F(T(p u ), p u ), F := sup F (T), (5) 

Tee 

F(T):= ( dwq{u)F{T{ Pul ), Pul ), F:=supF(T). (6) 
J Tee 

We will restrict attention to ensembles {p u } with a fixed degree of squeezing s, so 
that the benchmarks will be functions of s and we will occasionally write F(s), Fq(s) 
to emphasize this dependence. 

As usual, for the ideal scenario of pure states we use the fidelity 

F(T{p u ), p u ) = Tr[T(p w )p w ]. Clearly, in practice it is more realistic to assume 
that the initial pure states p^ undergo a noisy channel Af and become mixed before 
entering a quantum memory or a teleportation scheme. A possible option for F 
would then be the "Uhlmann fidelity" [53] between two mixed states, given by 



/(p 1; p 2 ) = Tr \J p\ /2 P2p\ /2 ■ Major drawbacks of this choice are that its non- 
linearity leads to a very involved theoretical optimization and that such a quantity is 
exceedingly difficult to measure in experiments (without invoking a plethora of extra 
assumptions) . We will present later on an analytical result adopting Uhlmann fidelity, 
but we will mainly follow a different route and use instead the overlap 

F(T( PuJ ), PuJ ) = Tr [T(N{p u )) Pu ]. (7) 

This quantity is easier to determine in experiments (p u may be interpreted as an 
observable) and by definition Fq and F are proper benchmarks in the sense that 
beating their values means outperforming any classical scheme. Strictly speaking, of 
course, this overlap measures how close the output of T is to the initial noiseless state 
rather than to the input. 

The derivation of the benchmarks simplifies considerably if the probability 
measure q is invariant with respect to a symmetry group G, i.e., if q(w) = q(g(uj)) for 
all g £ G and p g (u) = Ugp^U^ for some unitary representation U g of G. A standard 
argument [T3l [25] then implies that, w.l.o.g., the channels T in the optimization of 
Fq and F can be taken covariant with respect to G, in the sense that, for all density 
operators p and all g S G, one has 

U g T(p)Ul=T(U g pUl). (8) 
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For a compact symmetry group (with Haar measure dg) the argument is straight 
forward since we can replace every single T by a covariant counterpart 

f(p) := J dgUlT(U gP Ul)U g , (9) 

which performs at least as well as T if the chosen functional F is concave, as is the 
case for all the instances we discussed above. Moreover, if T is an element of a convex 
set closed under the group action, like the set of entanglement breaking channels, 
then so is T, If the orbit of G covers the entire parameter space then F(T(p u ), p u ) 
becomes state independent and we find Fq = F, i.e., the average case and the worst 
case performance become the same. 

The argument becomes more subtle if G is not compact, as in the case where no 
a priori information is given about the displacement £. In this case Fq is still well 
defined but already F has to be discussed more carefully, since it formally requires an 
average over a 'flat' distribution in phase space. Nevertheless, an analogous argument 
goes through and we can w.l.o.g. take T to be phase space covariant, i.e., for all £ and 
density operators p 

W 6 T{p)Wl = T(W 6 pW\). (10) 

Clearly, due to non-compactness one cannot make the averaging procedure explicit, 
but has rather to invoke an invariant mean whose existence is guaranteed by the axiom 
of choice. For a more formal discussion of these matters see [57] . 



4. Analytical benchmarks 

4-1. Classical benchmark for pure squeezed states under uniform rotations and 
displacements 

Let us now focus on the specific case of an initial single-mode squeezed state p s with 
CM 




which is uniformly displaced and rotated in phase space. The input ensemble is given 
by Pu = W^UepsUlW^ , where U$ = exp[i#n] denotes a phase space rotation (while 
W(: is a displacement operator, defined in Section [2]), and no a priori information is 
given about 8 6 [0, 2ir] and £ € R 2 , whereas the degree of squeezing s is fixed and 
known a priori. 

In order to determine the best possible fidelity achievable by a measure-and- 
prepare scheme on such an ensemble, we will first prove an upper bound by enlarging 
the set of allowed channels to the set T D £ of time-reversible channels. That is, 
to all the channels which remain valid quantum operations when concatenated with 
time-reversal (this includes, for instance, all measure-and-prepare schemes which are 
assisted by PPT bound-entanglement shared between sender and receiver). In a second 
step we will show that the obtained bound is tight, even within the set 5, as it turns 
out to be achievable by heterodyne measurement followed by preparation of a coherent 
state. 



Theorem 1 (Benchmark — pure states) Let T be either the set of all measure- 
and-prepare schemes or the larger set of time-reversible channels. Within these sets, 
the maximal achievable fidelity for the teleportation or storage of squeezed coherent 
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states with squeezing s and subject to uniformly distributed rotations and displacements 
in phase space is given by 



F(s) = F (s) = sup inf Tr [ PuJ T ( Plj )} 



1 



(11) 



Proof We first prove the upper bound. As already indicated above we can restrict 
ourselves to the set T of phase space covariant channels satisfying 

f(W 6 AW^) = W i f{A)wl, (12) 

for all Weyl operators W s , £ e M 2 and all A <E € 1 (L 2 (R)). This leads to F = F and 
lifts the need for the minimization over displacements HZ] • That is, we are left to 
determine 

UePsUlfiUePsUl) 



sup 



inf Tr 

?e [o,27r] 



(13) 



(in other words, the infimum is only taken over all squeezed coherent states centred 
at the origin, i.e., with vanishing first moments). 

Furthermore, time-reversible phase space covariant channels have a particular 
simple form in the Hcisenberg picture, where we will denote the map under 
consideration by T* . By Lemma |4] (see |Appcndix A[ ) they act on Weyl operators 



as 



where 



t*(w ( ) = 



t(0 = Tr 



tW, 



(14) 



(15) 



for some density matrix r. Conversely, every such r yields an admissible phase space 
covariant channel T . This allows us to recast the optimization over channels into one 
over density operators, as follows. 

Exploiting the Parseval relation ([3]) back and forth, we get for any density matrix 
p whose characteristic function is a centred Gaussian: 

PT(P)] = — / d 2 £Tr [pW s ] 2 Tr \tW. 



Tr 



2tt 

- i l d2 « Tr 



V2£ 



Tr 



tW 



V2£ 



(16) 



(where we took advantage of the fact that the characteristic function of centred 
Gaussian states is real and has a purely quadratic dependence on £). 

To compute Eq. (|13p it is now convenient to reintroduce an average over rotations 
(instead of an infimum, which gives the same value due to the optimization over T), 
and write Eq. (fT3"l) as: 



sup Tr 



1 

An 



2 K 



- u BPs uue r 



1 



271 



— / Ueps U<d9 



(17) 



This operator norm (largest eigenvalue) can be promptly determined since averaging 
over rotations just sets all off-diagonal elements in Fock basis to zero, so that one 
simply has to resort to the expression of a squeezed vacuum state in Fock basis, given 

by 

/Vi" y^Ql f 8-1 



l + s 



71=0 



2s + 2 



\2n) 
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As the largest diagonal entry is the vacuum component we finally obtain 

F(s)<Uo\p s \0) = 



2 1 lr " ' 1 + S ' 

where the r.h.s. is attained within the set of time-reversible channels. 

To show that this upper bound is actually tight for measure- and-prepare schemes 
as well, we will now explicitly point out that a specific, simple measure-and-prepare 
scheme attains the bound [28 . Consider a heterodyne measurement, i.e., a POVM 
{W^|0)(0| where the outcome £ is followed by the preparation of a coherent 
state Wf |0), so that the whole map H reads 

H{P) = j d 2 e (0| WlpWz |0) Wi |0) (0| W\ (18) 

This process acts on covariance matrices as 7 1— *■ 7 + 2 • 1 (while leaving first moments 
unaffected), so that one can easily determine the input-output fidelity for any state 
in our ensemble using again Parseval's relation ([3|) and solving the resulting Gaussian 
integrals [29] 

Tr \p u H{p w )] = 2 [bet ( 7ps + (j Ps + 2 ■ 1))] = (19) 
which is independent of £ and 9 and achieves the above bound. □ 

Note that the optimal strategy (heterodyning followed by generation of coherent 
states) does not depend on the degree of squeezing and is in fact the very same strategy 
which is optimal for a flat distribution of coherent states [5] (for coherent states, such 
optimality extends to rotationally invariant, Gaussian distributions of displacements 
as well [H])- The optimal fidelity as a function of the squeezing s is plotted in Fig.Q] 

This analytical result, though obtained for an ideal case, fully highlights the 
importance of radomising the input phase in practical instances. In fact, for a 'flat' 
distribution of displacements but fixed phase, the benchmark can be promptly inferred 
from the coherent states' case, and is just 1/2 > i/s/(l + s), regardless of the 
degree of squeezing s assumed. Randomising the input phase appears to be very 
helpful to provably reach the quantum regime in experiments (and, not surprisingly, 
the advantage granted by random rotations becomes more and more relevant with 
increasing squeezing): for pure states and degrees of squeezing within experimental 
reach (s < 10), the difference is around 0.2 and might easily turn out to be crucial for 
experimental success. As we will see later on (see Sec. l5.3p . and again not surprisingly, 
the same holds true for the randomisation of displacements. 



4-2. Classical benchmark for mixed squeezed states under uniform rotations and 
displacements 

Let us now consider the case of an input ensemble of mixed Gaussian states, derived 
from the initial pure states p u , lu = (s,£,0) by the application of a noisy Gaussian 
channel M, adding classical Gaussian noise with variance 77, and thus acting on 
covariance matrices as 7 1— > 7 + rjl while leaving first moments invariant [21] . Af 
can be understood as the application of random displacements to the input state, 
distributed according to a Gaussian with variance rj. We will keep referring to s as 
squeezing parameter (here, prior to the action of noise) and will refer to r\ as to the 
'additive noise'. 
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Figure 1. Fidelity benchmarks F(s,r]) for measure-and-prepare schemes on 
ensembles with squeezing s, flatly distributed displacements and random phase 
space orientations, as a function of the squeezing s (before additive noise is 
applied). The continuous (red) curve refers to pure states (rj = 0) while the 
dashed (blue) and dotted (green) curves refer to mixed states, for r\ = 0.5 and 
rj = 1 respectively. Notice that, as one should expect, the pure state case is an 
upper bound for mixed states for the figure of merit F. As stated in Theorem [3] 
the pure states case is also an upper bound for the Uhlmann fidelity benchmark 
of such mixed states. 



As pointed out earlier we will use the overlap with p u as figure of merit: 

Fo(s,r]) = sup inf Tr 
Tee C>0 




where, due to the symmetry (which is preserved by M), one has again Fq = F. For 
the same reason we can again restrict to phase space covariant channels and proceed 
along the lines of Theorem [TJ 



Theorem 2 (Benchmark — mixed states) Let T be either the set of all measure- 
and-prepare schemes or the set of time-reversible channels. Within these sets the 
maximal achievable overlap for the teleportation or storage of mixed squeezed coherent 
states with squeezing s and additive noise rj is given by 



F(s, V ) 



sup inf Tr 



T 



1 



2 +S 



-1/2 



(20) 



Proof As for Theorem [TJ we restrict to phase space covariant channels T 
parameterized by a density operator r and obtain in an analogous way 



Tr 



^(A.) 



Pa. 



1 

2^ 



d 2 £Tr 



tW. 



V2£ 



■fc« = ±T>[r//] 



(21) 
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where p' is now a centred Gaussian state with covariance matrix T = 7 p , +77/2. 
Once again, the supremum over r and thus over all time-reversible channels can be 
calculated by considering the state p' in the number basis, yielding 



sup— Tr 

T 4tt 



MUep'Ul 



1 

4~ 



(n|p'|n) 



:(0|p'|0) 



1 
2 



-1/2 



which is Eq. (|20p for time-reversible channels. Again, it can be shown by direct 
evaluation that the heterodyne strategy described by the map H of Eq. (fT8|) attains 
this bound, which completes the proof. □ 

Notice that the overlap achieved by the ideal quantum channel (identity) is 

-1/2 



(3 + l) (* 



Also, F(s, 0) = -F(s), as it should be since both the input 
ensemble and the figure of merit coincide with those of the preceding section in the 
noiseless case (77 = 0). Clearly, the optimal classical performance degrades with 
increasing noise (this is obviously the case for quantum strategies as well). Examples 
of such benchmarks for mixed states are displayed in Fig. [TJ 

For the same ensemble of mixed state, we shall also present a further analytical 
benchmark in the form of an upper bound, this time adopting the Uhlmann fidelity 
f(T(p),p) as figure of merit: 



Theorem 3 (Benchmark — mixed state fidelity) Let T be either the set of all 

measure- and-prepare schemes or the set of time-reversible channels. Within these sets 
the maximal achievable worst case ( or average ) Uhlmann fidelity for the teleportation 
or storage of mixed squeezed coherent states with squeezing s and arbitrary additive 
noise channel AT is bounded as follows 



supinf f(T{N{ Pu )),N{pJ)) 



< 



1 



Vry e R. 



(22) 



Proof Concavity of the fidelity allows us again to restrict to phase space covariant 
channels. All these channels commute with channels of the form Af*(W^) = 
W% exp[— 1 1£| \ 2 rj/4] (AT* standing for channel Af in Heisenberg picture). The result 
follows then from the contractivity of cp-maps with respect to the fidelity, together 
with the pure state result of Theorem [1] [Eq. fTTj) ]. □ 



5. Quantum benchmarks derived by Semidefinite-Programming 

5.1. Problem settings 

So far, we treated input ensembles of squeezed states which are displaced in phase 
space according to a 'flat' distribution. Needless to say, this is an idealization as one 
cannot implement an arbitrarily large displacement in practice. 

To be more realistic we have to treat an input ensemble of squeezed states whose 
first moments are essentially contained in a finite region of phase space, e.g., due to a 
sufficiently rapidly decaying probability distribution as in [5J [S] ■ 

This section will deal with this scenario by resorting to numerical means, as 
a purely analytical treatment appears to be far too involved for finite, non-flat 
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distributions of displacements (where the restrictions to phase space covariant channels 
is no longer optimal). We will show that the problem of computing benchmarks of 
the desired kind can be cast into a semidefinite program (SDP). As such, it comes 
with a guarantee of computing the correct value, as SDPs come in pairs of a primal 
and a dual problem which yield converging upper and lower bounds to the sought 
solution [30j . As the original SDP is in an infinite dimensional space, truncation will 
be necessary and will induce errors. However, we will provide a rigorous bound to 
the truncation errors so that the finally derived benchmarks are reliable and rigorous, 
constituting upper bounds to the actual optimal classical figures of merit. 
The figure of merit we use is again the average overlap 

F = sup / dijj q , (w)Tr 
TesJ 

where Af is a noisy channel that describes the noise suffered by PuJ before entering the 
storage or teleportation device, which for our purposes should be 

(i) a channel which allows for the computation of the matrix elements in the Fock 
basis (k\N{pJ)\l), e.g., a Gaussian channel [H] 

(ii) rotationally covariant, i.e., J\f (Ugpuf) — UeN (p) Ug for all p, 9. 

Fortunately, channels representing attenuation, amplification, and thermal noise are 
all of this type, the attenuation channel being probably the most relevant here, as it 
models losses of photons/excitations which are the dominant decoherence process in 
the practical realisations we have in mind, involving traveling waves of light. Thus we 
will adopt it in the following, denoting it by A/a. The parameter A € [0, 1] represents 
the transmitivity (in intensity) such that the channel acts on, respectively, covariant 
matrix 7 P and vector of first moments d p as: j p i— > + (1 — A)l and d p i— > ^/\d p . 
Another crucial assumption for our method is that the input distribution q is uniform 
in the angle 9. Moreover, we will consider ensembles with constant degree of squeezing 
s (although this is not necessary for the method). Hence q can be considered a 
probability distribution which depends only on the displacement £ and therefore 

F = sup / / g(£)Tr 

Tee J0e[o,27r] J^s. 2 

where p s ^ is a pure squeezed states with degree of squeezing s and displacement £. 
Note that F in this way becomes a functional of A, s and of the distribution q. 



Tbf( Pu y\ Pl 



T(M x (UePs,tUt)\(u gPa ,tUn —d^ (23) 



5.2. Reduction to a finite dimensional semidefinite program 

In this subsection, we will see how we can reduce the quantum benchmark F in Eq. (|23f 
to a finite dimensional SDP. We first show that Eq. (|23|) can be reduced to an infinite 
dimensional SDP problem. To this end, we use a simple correspondence between the 
set of all entanglement breaking channels and the set of bipartite separable positive 
operators [31 which is nothing but the Choi-Jamiolkowski isomorphism |32j . albeit 



for an infinite dimensional system (see Appendix B for a proof) 



Theorem 4 (Choi-Jamiolkowski) Suppose 03 (W) is the space of all bounded 
operators and <£i(H) is the space of all trace class operators on a separable Hilbert 
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space TL. Then, for all entanglement breaking channels ^> onTL, there exist a unique 
separable positive bounded operator onTL®TL such that TrB(fl(^)) = 1 A and 

Tr(BV(A)) = Tr(n(^)A <g) B) (24) 

for all A € <Li{Ti) and B G 93(7i). Conversely, for a separable positive bounded 
operator Q, onJi^Ji satisfying Trs(^l) = 1a, there exists a unique channel ^(Q) 
such that it is entanglement breaking and satisfies Eq. 



By means of the above theorem, we immediately derive an upper bound to the 
quantity F of Eq. (|2"3|) in an infinite dimensional SDP form by enlarging the set of 
positive separable operators (denoted by Sep) to the set of positive operators with 
positive partial transpose (PPT) fi r > 0: 



def 

V = 



F = sup |Tr(fh7) ft € Sep, Tr B ft = 1 A \ 

< sup {Tr(«77) > 0,n r > 0,Tr B n = 1 A ) . (25) 
In the above equations 77 is a state on Tt <g) Tt defined by 

f q(0 UeUxiPsM ® (UePstUt) f- d£, (26) 

6»G[0,27r] J^m 2 v / Z7T 

where we exploited the rotational covariance of J\f\ . 

Here, we should remark that, if A = 1 and q(£) = — exp[— a\ |£| | 2 ] for a > 0, the 
inequality in Eq. (|25p is actually an equality. This is because, under such assumptions, 
we can choose an optimal f2 to be Gaussian [5] and the PPT condition f2 r > is 
necessary and sufficient for a two-mode Gaussian state fi to be separable [33J 131] • 

In the remaining part of this subsection we will transform the above infinite 
dimensional SDP into a finite dimensional SDP. We will need the following statement: 



Lemma 1 (Operator norm bound) If a positive separable operator il £ <8(7i<g)7i) 
satisfies Tr B fi = t A , then, £1 also satisfies \\fl\\oo < 1j where || ■ is the operator 
norm. 

Proof From the proof of Theorem 2] can be written as fl(^f) = J x M(dx) ® a(x) 
by using a POVM {M(dx)} and a set of states {a(x)}. Then, for any normalized state 
|$) on Tt <g> Tt, we can bound ($| Q |$) as follows: 

= / Ti(M(dx) ® o~(x) |$) ($1) 
Jx 

< [ Tr(M{dx)p A ) = l, 
Jx 

where p A is defined as Try3(|$) ($|) and we used that a(x) < 1 together with 
M(X) = 1. Therefore, we have \\Cl\\oo = sup^, ($| Q |$> < 1. □ 

Notice that, even if HOUoo < 1 was shown for the corresponding operator of an 
entanglement breaking channel, this fact may not be true for other channels (in general 

we can only show that the geometric measure [35] G(fi(T)) = svLp aeSep Tr (O(T)ct) < 
1). 
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Before we derive a finite SDP problem from Eg. ([25)1 . we observe the following 
crucial fact: By means of the action of the group integral over {Ug £§) t/e}ee[o.27r]i V is 
block diagonalized as 



V 



c=0 



Qc 



9(0 MO^e) ® d£ 



(27) 



where Q c d = £ fe+z=c |fc) 

dcf 



® \l) (l\ is a finite dimensional projector. Thus, defining 

Pc == X)i=o Qi> we obtain rj — P c r)P c = YliL c +i QiVQi > for all c. We arrive at 
the finite SDP upper bound now as follows: Suppose a separable positive f2 satisfies 
Treft = 1 A - Then, 

Tr(Or?) = Tr(r»F c r;P c ) + Tr{Q{n - ^P c )) 

< Tr(nP cV P c ) + HfilU Tr (77 - P c ryP c ) 

< Tr(f7P c ?7P c ) + Tr(ry - P c r]P c ) 

= Tr(f>P c 7/P c ) - Tr(P c ryP c ) + 1, (28) 

where we used the positivity of r\ — P c r]P c and Lemma [T] 

Finally, by means of Eq. ((28|) . we can derive an upper bound to the first line of 
Eq.([25|). Denoting with R c — (J2i=o I*) ® (J2i=o N) tne projection onto the 
subspace with a photon number smaller than c in each mode, with support suppP c , 
we obtain 



F < 



sup 

fie as (H) 



jTr(fip?7P c 



il E Sep, Tr B n = 1, 



< 



ne<B(H) 
sup 

SlG23(supp_R c ) 



< 



} + 1 - Tr(P cV P c ) 
sup \Tr(R c rtR c P c rjP c ) Q € Sep, Tr B fi = lA + 1 - Tr(P c rjP c ) 
^Tv(nP c riP c ) n € Sep, Tr B fi < 1a} + 1 - Tr(P c ?7P c ) 
{Tr(fLPriP c ) n > 0, fl r > 0, Tr B < I A + l - Tv(P c r]P c 



sup 



(29) 



The above upper bound now only involves a finite dimensional SDP, as the original 
infinite dimensional fl has been replaced by the finite dimensional R c flR c , with the 
same type of constraints (note that R c was introduced to preserve the separability of 
the operator). Similarly, the term Tr(P c ?7P c ) is just a trace of a finite dimensional 
matrix, so that we can numerically compute every term of the last formula of 
Eq. (|2"5|) for as large a truncation parameter c as our computer allows. Thus, Eq. (|2l))) 
enables us to efficiently calculate an upper bound for the quantum benchmark F for 
any probability density q(£), that is, ultimately, for any rotationally-invariant input 
ensemble. Moreover, since in the limit of large c we obtain the second formula of 
Eq. (l25j) . we can expect that this upper bound reaches the exact value for this bound 
for sufficiently large c. Finally, we rephrase Eq. (|29[) in the form of a theorem: 

Theorem 5 (Benchmark — SDP) With the above definitions, for any probability 
density g(£) on R 2 and rotationally covariant noise channel Af, we have 



F < sup 

rie<8(suppK c ) 



{Tr(OP c 7?P c ) n > 0, il r > 0, Tr B CI < 1 A | + 1 - Tr(P c ?7P c ) 



(30) 
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Figure 2. A plot of the known optimal classical bound F ■■ 



2a + l 
2{a + l) 



the upper 



bound Fi n fi n ite derived by our numerical calculation with c = 35, and the result 
of the finite SDP 

F finite (without corrected truncation error) for an ensemble of 
coherent states, no noise (A = 1), and q(£) = — exp[— a||£|| 2 ]. 



5. 3. Results of numerical calculations 

In this subsection, we present results of numerical calculations which were produced 
using inequality (|30[> . In order to reduce the memory requirements, we further exploit 



the block diagonal structure of rj in Eq. ([27f and impose J2j=o Qj^Qj — ^ f° r the 
implementation. The evaluation of the bound then splits into two terms: 

F finite d = sup {Tr(fiP c 7?P c ) > 0, n r > 0, Tr B ft < 1 A ) , 
ne<8(su PP fl: e ) 1 J 

terror — 1 Tr(_P c ?]_f , c ) . 



Denoting the r.h.s. of (piO)) by F in fi n it e we have 

F in finite — F finite ~t~ terror- (^1) 

Our computation proceeds along the following steps: 

(i) Compute matrix elements of J\f\{p s ,{) \l\) and {k2\p s ,(, ^2) for all ki, fe, '1, 
I2 satisfying ki + fc 2 < c and l\ + I2 < c; c is a fixed maximum photon number 
practically bounded by the computer memory. For the calculation, we used an 
analytical formula for matrix elements of single-mode Gaussian states as a finite 
sum over Hermite polynomials (see Eq.(4.10) of 36J ) . 

(ii) Evaluate the following integral over £ 

/ q{0{ki\mp S ,d\h)(k2\p*M 
by Quasi-Monte Carlo with the Halton sequence [57] . 
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Figure 3. Results of numerical calculations of the quantum benchmark Fi n fi n i te 
for an ensemble of noisy squeezed states displaced according to a Gaussian 



probability distribution q(^) = — exp[— a\\ 
squeezing s = 8, and losses (1 — A) = 0.4, 0.2, 0. 



Levels used in Fock space c = 30, 



(iii) From Ea.([2?|). we derive P c r]P c . 

(iv) Finally, we solve the SDP leading to -Fin/mite- 

Here, all the above numerical calculations were implemented in Matlab, and the SDP 
is solved using the Matlab toolbox SeDuMi version 1.1 [35]. In the above numerical 
calculation, we do not use any approximation except for the Quasi-Monte Carlo 
integral. 

As a first application, in order to provide the reader with a convincing test for our 
numerical calculations, we show that our method reproduces the result of the optimal 
quantum benchmark for displaced coherent states derived in [S]. Suppose A = s = 1 
and the input ensemble is distributed according to = — exp[— a||£|| 2 ]. That is 
we deal with an ensemble of coherent states whose centers are distributed according 
to a Gaussian distribution with variance l/(2a). Hence, using the invariance of the 
ensemble under rotations we obtain 

V= I -exp[-a||£|| 2 ] p u ®p u d£. (32) 
Jim 2 71 

The quantum benchmark F in this case was shown [9] to be (note that our definition 
of the parameter a is different by the factor 2 from the definition in [5]) 

— 2a + 1 

2(a + 1) 

So we can compare our upper bound F \ n finite derived by SDP with the optimal bound, 
which is shown in Figure [2] for the maximum photon number c = 35. As expected the 
numerics satisfies Fi n fi n it e > F > F finite but we observe in Fig. [2] that as long as 
a > 0.2, the values of both Fi n fi n i te and -F/mite are essentially indistinguishable from 
the optimal bound F while for a < 0.2 we find noticeable differences. In other words, 



Quantum benchmarks for teleportation and storage of squeezed states 15 
Quantum benchmarks randomly rotated squeezed vacua 



o 



o 

o 

o 

o 

o 

o 

o- 

■ a 

123456789 10 

so 

Figure 4. The results for the quantum benchmark F in j inite for an ensemble 
of randomly rotated squeezed vacuum states without displacement, i.e., q(^) = 
5(| |£| | 2 ), c = 30 and no noise (A = 1). 

in Eq. (|31|) the term e err or, which originates from the truncation of the dimension, 
becomes dominant in this region. As a result, the minimum value of fin/mite is 
around a = 0.15, i.e., -Fm/imte is not monotonically decreasing with respect to a. 

From this result, we may expect that our upper bounds are also almost optimal 
in other situations (different q, s,Af) as long as e error is small with respect to -F/mite- 

Next, we consider the case of an ensemble of noisy squeezed states with a fixed 
squeezing parameter s and again a Gaussian prior distribution for £. We show the 
results of our calculation in Figure [31 Here, we chose c = 30, s — 8, and plot for 
A = 0.6,0.8,1. In this figure, we can observe that, as expected, Fi n fi n i te decreases 
with decreasing A. Due to the increasing contribution of terror for decreasing a, the 
best (lowest) values of our benchmarks are achieved, under realistic losses, for a ~ 0.15. 
This value, though maybe non-optimal, is suitable for comparisons with realistic 
experimental situations. A simple, preliminary analysis of possible experimental noise 
conditions indicates that the values of the benchmark for such parameters could be 
beaten by forthcoming experiments aimed at the teleportation or storage of squeezed 
states. Let us mention that single states have already been teleported or stored, 
both with light modes [TSJ [Hj and with atomic memories (TTJ, [TS] . The experimental 
demonstration of the transmission of a full ensembles of squeezed states is yet to come, 
but techniques are ripe for it to be realized in both settings: our method is ready for 
the analysis of such developments by direct comparison. 

Finally, we consider the simple case of an ensemble of randomly rotated squeezed 
vacuum states without displacement in phase space. That is, we choose g(£) = 
(S(| |£| | 2 ), shown in Figured] Here, we chose c = 30 and plot for 2 < s < 10. Note 
that for s = 1 there is only one state in the ensemble, the vacuum state, so that 
F = 1 in this case. Fig. [3] clearly shows that, if the displacements in phase space are 
not randomised, the classical benchmarks we derived increase significantly and and 
beating them in current experiments would be a daunting challenge: distributing the 
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displacements of the input ensemble is thus definitely, at present, a technical necessity 
in order to achieve a certifiable quantum performance [39] . 

6. Summary and Conclusions 

In this work we have derived upper bounds on the performance of classical 'measure 
and prepare' protocols for the quantum storage and teleportation of ensembles of 
squeezed states of continuous variable systems. These bounds may be employed as 
benchmarks to discern truly quantum mechanical performances in experiments aimed 
at the realization of such processes. Motivated by currently available experimental 
capabilities and concrete set-ups, we have considered ensembles comprised of states 
with fixed squeezing, but allowed for the possibility of random rotations and 
translations of those states in phase space. 

Fully analytical benchmarks have been obtained for pure squeezed states as well 
as mixed states obtained from the application of additive thermal noise, for the case of 
completely unknown rotations and translations in phase space. In these ideal cases, the 
benchmarks decrease monotonically with increasing degree of squeezing s, vanishing 
for diverging s: in this sense, increasing the squeezing reveals the infinite dimensional 
character of the Hilbert space. Furthermore, we have presented a numerical technique 
that first reduces the problem to a finite dimensional setting with rigorous error 
bounds, and then allows for the numerical solution of the remaining finite dimensional 
problem using semi-definite programming. The so obtained bounds are rigorous and 
this approach may be applied to more general settings, even beyond the Gaussian 
regime, for variable squeezing and arbitrarily distributed displacements, with the only 
proviso that the distribution over rotations should be still rotationally symmetric in 
phase space, to maintain covariance under rotations (which was essential in deriving 
rigorous bounds for the truncated problem). Let us also emphasise that our numerical 
strategy can be reliably applied to finite sets of states as well (rather than to continuous 
distributions), which are what, strictly speaking, is sampled in actual experiments. 

Our numerical results strongly indicate that allowing for randomised displace- 
ments is crucial to outperform optimal classical strategies with current quantum tech- 
nologies and that, on the other hand, if such a randomisation is allowed, realistic 
setups might be able to enter the quantum regime with squeezed states. Likewise, 
randomised phases, other than providing a much more "appealing" evidence for the 
teleportation or storage of the states (since, in such a case, states with a varying struc- 
ture of second moments would be transmitted), seems to be needed for the benchmarks 
to be beaten. 

Overall, the results presented here will provide a useful resource for assessing 
the quantum mechanical character of experiments aimed at demonstrating quantum 
storage and transmission of squeezed states. In a future work, these methods will be 
applied to concrete experiments. 

Note added. We acknowledge that, very recently (after the completion of this study 
and during the writing of the present paper), another work appeared (Ref. [12] ) where 
benchmarks based on the Uhlmann fidelity for rotationally covariant ensembles were 
studied, and the reduction of the benchmark estimation to an SDP was independently 
derived. 
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Appendix A. Phase space covariant channels 

In this appendix, we collect useful results about 'phase space covariant' channels, i.e. 
channels covariant under the action of the Weyl-Heisenberg group of displacement 
operators G R 2 ™}. For a phase space covariant channel one has 



T(Wj<)W e ) = W^T(0)W^ , VO G Q5(Z, 2 (R)® n ) and V£ G 



1)2 n 



[03(^2 (R)®") being the set of bounded linear operators on n copies of the bosonic 
Hilbert space]. These results were already used (though not explicitly detailed) in 
[2"5I 121)] . Central to this analysis is the characterization of the class of "linear bosonic 
channels" in terms of their action on Weyl operators, given in [JD], which we report 
here without proof in the form of the following lemma. Recall that T* stands for the 
operation T in the Heisenberg picture. 

Lemma 2 (Linear bosonic channels) A map T*(W^) := f{t,)Wx(, is a quantum 
channel if and only if / is the quantum characteristic function with respect to a 
modified symplectic form a := a — X T aX. That is, / : R 2 ™ — + C has to be continuous, 
/(0) = 1 and every matrix with entries 

M H = /(^ fe )-e W )exp[i^.^«] (A.l) 

has to be positive semi- definite for all £( fe ),£W G R 2 ™. 

Note that condition (|A. 1[) for a = a (a — 0) is equivalent to / being a quantum 
(classical) characteristic function. 

The following characterization of phase space covariance ensues: 

Lemma 3 (Phase space covariant channels) Every phase space covariant chan- 
nel is uniquely characterized by a classical characteristic function / and acts in the 
Heisenberg picture as 

T*(Ws)=f(0W 6 . (A.2) 

Proof Consider the action of T* on a Weyl operator W%. Exploiting the Weyl 
relations, phase space covariance (with respect to any W v ) leads to 

W v T*(Wt) = e^" 7 " T*(W^)W n . (A.3) 

It is straightforward to check that this entails 

[T*(Wt)Wl,W v ]=0, V^GR 2 ". (A.4) 

The irreducibility of the Weyl system then implies T*(W^)W} oc 1. Denoting the 
proportionality constant by /(£), we have a map of the form in Lemma[2]with X = 1. 
Hence, / has to be a classical characteristic function. □ 
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Let us denote by d the time reversal (or matrix transposition) operator. Every 
entanglement breaking channel, i.e., 'measure-and-prepare scheme' (see Ref. |23j). is 
such that T o <d is completely positive (that is, the Choi matrix - or the Jamiolkowski 
state - has positive partial transpose). Moreover, time reversal acts very simply in 
phase space by flipping the sign of one of the two canonical quadratures (see, e.g., 
[33]): #(Wf) = Wz£, where Z = ©™ =1 cr 2 (a z being the Pauli z matrix). We will now 
apply this additional constraint to phase space covariant channels, in order to achieve 
a stronger characterization: 

Lemma 4 (Time-reversible channels) A phase space covariant channel T is such 
that T o is completely positive iff it has the form 

T*(W £ )=/(e/V2)W £ , (A.5) 

where / is any quantum characteristic function, i.e., there is a density operator 
r e £i(i 2 (R) c5 ™) such that /(£) = Tr [tW{\. 

Proof Using Lemma [3] we obtain 

$oT*{W t ) = f(0#(We) = f&Wzt , (A.6) 

where Z reverses the momenta so that Z T aZ = — a. Following Lemma [3 this map is 
completely positive if and only if positivity of Eq. (|A.1[) holds for a — 2a. The result 
follows then by substituting £ — > £/v2- D 

Appendix B. Proof of Theorem [4] 

Proof Suppose W is a entanglement breaking channel and can be written [23j [22] as 
*G°) = fx Tr(M(dx)p)o-(x) for all p £ £1 (ft), where a(x) E £i(H) are states and M 
is a POVM, i.e., M(X) > for all Borel subsets A of a complete separable metric 
space X for which M (X) = 1. 

Then, = f x M(dx)®a(x) satisfies Tr s (fi(*)) = 1 A and Eq. (gj) . Moreover, 

suppose there exists a 0'(*) satisfying Tr B (Q'(^)) = 1 A and Eq.([2I]). Then, 
since Tr(fi'(*) |t) (j\ ® |fc) = Trflfc) (i| (j|)) = Tr(fi(*) |t) (j\ ® (Z|) for 
a orthonormal basis of 7i, we immediately have il' = £7. 

Conversely, suppose there exists a separable positive bounded operator f2 on 
H®H satisfying Tr B {Q) — 1a and can be written as = / |$)($| <8> I^K^I p.(d$dtp) 
with measure /i. We define a channel \& as p i— > J Tr(M(dy;)|o)|y)((^|, where 

M(#) = f J $ |$)($|m(c^#0- Then, since Tr B (fi) = 1 A , M is a POVM; that is, 
'J is an entanglement breaking channel. Evidently, \& satisfies Eq. ([2"4"]) . Moreover, 
suppose an entanglement breaking channel ^' also satisfies Eq. ([24]) . Then, since 
Tr^l-'LA)) = Tr(fi(#)A<g) B) = Tr(£#(A)) for all £ € «8(ft), tf'(A) = for all 

A e £i(7i)- □ 
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